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Abstract 

The plane case of central configurations with four different masses 
is analyzed theoretically and is computed numerically. We follow 
Dziobek's approach to four body central configurations with a direct 
implicit method of our own in which the fundamental quantities are 
the quotient of the directed area divided by the corresponding mass 
and a new simple numerical algorithm is developed to construct gen- 
eral four body central configurations. This tool is applied to obtain 
new properties of the symmetric and non-symmetric central configu- 
rations. The explicit continuous connection between three body and 
four body central configurations where one of the four masses ap- 
proaches zero is clarified. Some cases of coorbital 1+3 problems are 
also considered. 

Keywords: Newtonian Four-Body Problem. Central configura- 
tions. 

PACS 45.50.Pk Celestial mechanics 95.10.Ce Celestial mechanics 
(including n-body problems) 
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1 Introduction 



In this paper we study the central configurations that appear as partic- 
ular solutions of the diflFerential equations of motion of the Newtonian 
4-Body Problem, namely [1] 

- ^ °"""'r^ = 2. 3. 4.) (1) 

where r^, (j = 1, 2, 3, 4) denote the position vectors of the particles in 
three dimensional space; rjj = rji are the distances between particles 
i and j; rrij, {j = 1, 2, 3,4, ) denote the constant masses, and G is the 
(measured) gravitational constant. The masses of the four bodies mi, 
1712, 1^3 and 7714 are generally different. The right hand side of this 
equation is the gravitational force which is derived from the poten- 
tial energy which is the Newtonian function of the distances between 
particles r^j. 

For any configuration to be central, vector forces must be in the 
direction of the position vector so that the right hand side of equation 
^ divided by the mass rrii equals a constant by the position vector 

The relevance of central configurations is remarked in multiple sites 
in the literature. We appreciate specially the beautiful presentation 
of the general theme by D. Saari [2j where the fundamental references 
are found. Central configurations are the only known cases where the 
differential equations of the N-Body Problem are integrable. They are 
relevant for the simultaneous collision of the N bodies, and in the limit 
that the relative positions of N colliding particles increase without 
bound. Topological reasons ask for their properties to be discovered. 
Even the number of essentially different central configurations is a 
mathematical challenge for this century. An additional new reason is 
their usefulness to demonstrate the existence of non integrability of 
the N-Body differential equations [3]. The original sources may be 
found in these two references. 

The study of central configurations of the N-Body Problem is a 
subject that expands in many important directions. Since central 
configurations of the Three-Body Problem are better known, we tackle 
the Four-Body Problem. 
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The Four-Body Problem in the plane was already considered by 
some researchers in the Nineteenth Century. Again we cite Albouy et 
al. [3] to find some of them. In 1900 Dziobek [4j published a seminal 
work on central configurations with relative distances as coordinates. 
After him remarkable papers were published, among which we select 
the MacMillan and Barker study [S], with different coordinates and 
similar results. Later Schmidt [6] reviews and extends these results 
for five particles in three dimensions. Moeckel [7] proves the finiteness 
of the number of Dziobek's configurations. Leandro [8] refines the 
finiteness problem and studies important bifurcations. 

Although the ideal for physicists and mathematicians is the general 
case of four different masses, the difficulties in simplifying it have 
lead to a prudent research based on many particular cases with some 
symmetry. The case of four equal masses was analyzed by Albouy [9]. 
When three of the masses are equal, but different from the fourth mass, 
a research by Bernat, Llibre and Perez-Chavela pj] settles interesting 
theorems on the number of solutions and bifurcation values in the case 
of kite symmetry. A reciprocal result to this is found by Long and Sun 
[llj . Properties of symmetric cases are also obtained by Albouy and 
coworkers [3]. 

Our main contribution to the Central Configurations literature is 
to set forth the four weighted directed areas as known parameters. 
Although the approach of not choosing the masses as fundamental 
data is against tradition, it is justified by the simplicity of the resulting 
algorithm that produces as output values for the six distances and the 
four masses in a compatible central configuration. 

This approach leads us to one non-trivial algebraic equation in a 
single unknown, to be compared with the many unknown variables 
approach used by others. 



The directed vectors representing twice the area of each of the four 
faces of a tetrahedron with vertices ri, r2, r^, r^, are: 



2 The plane problem 



Si 

S2 

S3 
S4 



r2 A ri + rs A r2 + ri A 



r4 A r2 + r2 A ra + rs A r4 
ri A r4 + r4 A ra + rs A ri 



ri A r2 + r2 A r4 + r4 A ri 



(3) 
(4) 
(5) 
(6) 
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The vectorial sum of these four vectors is the zero vector 



Si + S2 + S3 + S4 = 



(7) 



The case where the four particles at the vertices are in a constant 
plane is an important and old subject. The four vector directed areas 
are all parallel, the direction of which may be taken as the unit vector 
k. 

Sj = kSj , (8) 

where the Sj are twice the directed areas, positive if Sj is parallel to 
k, and negative if the direction of Sj is opposite to k. 

We choose the third component of the cartesian coordinates of the 
four particles as zero and form the matrix of cartesian coordinates 



X4 ^ 

y 

The four directed areas are written in terms of these coordinates 



Vi 



X2 
2/2 





3^3 





as 



S^ = 



53 



1 


1 


1 






X4 


y2 


ys 


2/4 


1 


1 


1 


Xl 




X4 


Vi 


2/2 


2/4 



^2 = 



1 


1 


1 


Xl 


X4 


X3 


2/1 


2/4 


2/3 


1 


1 


1 


Xl 


X3 


X2 


2/1 


2/3 


2/2 



(9) 



that are the four signed 3x3 minors formed from the matrix 

/ 1 1 1 1 \ 

Xl X2 X3 X4 (10) 

»^ 2/1 2/2 2/3 2/4 / 

Addition to the previous matrix of a row equal to any of its three 
rows produces a square matrix with zero determinant, implying that 
the necessary and sufficient conditions to have a flat tetrahedron are 



4 

1=1 



0. 



and 



i=l 



Xi 



0, 



4 

E 

i=l 



SiVi = . 



(11) 



(12) 
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The two last equations may be grouped in the zero vector condition 
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J2Siri = 0. (13) 



i=l 

From these properties we deduce a necessary condition for flat 
solutions in terms of distances. 

= E E ^S^Sj = E E(r? - 2r. • r,- + rj)S.S, = 

1=1 j=l i=l j=l 

2(X: Si)Ct - 2(X: S,r,) ■ Ct Sj^j) ■ (14) 

i=l j=l 1=1 j=l 

This is zero as a consequence of the previous conditions for a flat 
solution. Indeed the first term is zero because the sum of the directed 
areas is zero, while the second term is zero by equation (jl3p . Note that 
equations (fTTl - [T4l) are purely geometrical, independent of the origin 
of coordinates and independent of the masses. Beside the condition of 
the sum of directed areas equal to zero, other authors take into account 
the Cayley-Menger determinant, but we do not. Constraint ()14p will 
be seen to be very important in the theory of central configurations 
in section 5. Equations (jlip and (jlSp were considered by Dziobec [1] 
and are currently connected to affine geometry and the plane four- 
body problem (see Albouy et al. [3J and references to Albouy therein; 
note that equation (jl4p also follows trivially from equations used by 
Albouy who gives it [12j in explicit form.) 

Other important remark is that although other authors generalize 
results for n particles in n — 2 dimensions, and other forms of the 
potential energy, we do not follow that path. 

Implicit in what follows is Heron's formula relating the absolute 
value of the area of a triangle to the square root of a function of the 
three sides, which can be expressed, for example, as 



/ -1 1 1 \ / rfg 



1 -1 
1 1 



rh ■ (15) 



r 



3 Central configurations 

In this section we begin with Laplace's approach (see ref. [I3]). Cen- 
tral configurations are defined by the condition that the force divided 
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by the corresponding mass is proportional to the position vector: 



Br, = Y^'W_J^ (for i = 1, 2, 3, 4,) (16) 

where -B is a constant, the same for all the particles. Equations p6p 
imply that the origin of coordinates is at the center of mass. 

Taking the external product of equations i and k with (rj — rfc)A 
and equating, we find the six conditions 




in terms of the weighted vector areas Aj = Sj/mj 



Ai = 


1 r 

r4 A r2 + r2 A rs + rs A r4 

mi 


(23) 


As = 


[ri A r4 + r4 A rg + rs A ri] 

m2 

1 r 

ri A r2 + r2 A r4 + r4 A ri 

1713 


(24) 


As = 


(25) 


A4 = 


1 r 

r2 A ri + rs A r2 + ri A rs . 


(26) 



These equations are attributed to E. Laura and H. Andoyer (see ref 
i). 

In the non-planar case, obviously no pair of these four vectors 
are parallel for arbitrary masses and positions, and the unique three 
dimensional central configuration satisfying equations (I17ti22p requires 
all six distances to be the same, so that the masses lie at the vertices 
of an equilateral tetrahedron. This central configuration has the only 
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solution [13] of straight paths towards or away from the center of mass. 
In the rest of the paper we shah consider plane non-collinear central 
configurations of four masses. 

The plane solutions with zero enclosed volume but finite area are 
obtained by taking into account that the four vectors (|23ll26p are now 
parallel. They are written in terms of the weighted directed areas Aj 
defined as 

A,=Ajk, (27) 

with 

A, = ^, (28) 

where Sj is the value in vector Sj multiplying the vector k, equal to 
the directed area of the triangle with the three particles different from 
j at its vertices. 

Suppression of vector k from equations (117ti22p leads us to the 
homogeneous system 








^4 







Ai 


-A, \ 


( ^/rh 


\ 




( ° \ 




-Ai 





Ai 


-A2 





A2 


1/^31 











Ai 


-Ai 





A3 


-A3 





l/^?2 














A2 







A3 


-A2 


1/^41 











-Ai 





^3 


-^3 





Ai 


1/^2 









V 


Ai 


-A2 





A2 


-Ai 


1 


V l/r|3 


) 







(29) 



The matrix notation is ours. Since this equation may be considered a 
consequence of Dziobek's equations, they are frequently attributed to 
Dziobek. but they were obtained previously, in a similar form to that 
given here, by E. J. Routh (see [3j.) 

Since the coefficient matrix of this system is antisymmetric, a non- 
trivial solution exists that is the linear combination of its two eigen- 
vectors with eigenvalue zero, namely 

(111111) and {A2A3 A3A1 A1A2 AiAi A4A2 A4A3) . 

The first eigenvalue is evident considering the structure of the rows of 
the matrix and the equilateral solution. The existence of the second 
eigenvalue is secured because the matrix is antisymmetric of order six. 
This fact is the consequence that a real antisymmetric matrix becomes 
hermitian by multiplying its elements by the imaginary unit i. The 
eigenvalues of a hermitian matrix are real. Therefore, the eigenvalues 
of a real antisymmetric matrix are pure imaginary numbers forming 
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complex conjugate pairs or zero. If the order of the square matrix 
is an odd number we have at least one zero eigenvalue. If the order 
of the antisymmetric matrix is an even number, the zero eigenvalues 
appear in pairs. 

Hence the solution of (j29p is given in terms of two parameters A 
and a 



This equation was also obtained by Dziobek [3] from a different ap- 
proach to this problem, considering the critical points of the potential 
energy, constrained to a constant inertia moment and zero volume 
(this is the perspective of ref. [6]). One way to take into account the 
zero volume condition is by taking the derivative with respect to rfj of 
the Cayley-Menger determinant (which is proportional to the square 
of the volume of the tetrahedron). These derivatives are proportional 
to the product SiSj if the plane restrictions are taken in account. An 
obvious difficulty of that deduction is that since the non trivial entries 
of the Cayley Menger determinant are just the squares rfj, to obtain 
the square roots in the expression of the directed areas Si as functions 
of the rfj, equations (fT5|) . one needs to use explicitly restrictions ([TT]) 
and ()13p connecting these areas. Actually Dziobek's proof remained 
unpublished until the end of the Twentieth Century (see Moeckel [7J ) . 
Other proof using convex theory has been developed by Albouy (see [3] 
and references therein). Although D. Saari follows Dziobek's approach 
of critical points, he does not recommend [2] the use of Cayley-Menger 
determinant but instead he is in favor of geometric arguments. Note 
that our proof is quite elementary using simple properties of linear al- 
gebra and we follow Saari's recommendation of geometric arguments 
[2], using Dziobek's geometric restrictions of the previous section (jllh 
and p^ . also made explicit in Albouy et al. [3j. 

Now we obtain a well known expression for the parameter a, which 
starts by writing (j30p in the form 



Multiply both sides by and sum over j and k. Using the geometric 
plane constraint (jl4p we find 




a + XAjAk . 



(30) 



(31) 



X; mjnik/rjk 
j>k 



(32) 



a = 



J2 nijmkr'^-^ 



j>k 
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Figure 1: Two different arrangements of four particles with different masses 
following elliptic trajectories in a central concave configuration. The common 
focus of the four ellipses is at the right hand side focus of the small ellipse. 
The eccentricity of the four ellipses is e = 0.72. Ai = 5, A2 = 6, = 4, 
A4 = -8. 

expression that is positive definite. From either (j32p or (j30p it is seen 
that a has dimensions of the reciprocal of a cubed length. 

Further, substitution of (j30p in the equation that defines central 
configurations (jl6p . we verify that it holds with the constant B given 

by 

4 

B = -aJ2m^. (33) 

This just requires that the origin of coordinates is located at the center 
of mass and the plane constraints (fTTl) and (fT3]l . 

On the other hand, we have that A is negative. For the proof, we 
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1 



Figure 2: Two different arrangements of four particles with different masses 
following elliptic trajectories in a central convex configuration. The isolated 
point is the center of mass and the common focus of the four ellipses. The 
eccentricity of the four ellipses is e = 0.72. Ai = 15, A2 = —6, = 3, 
Ai = -4. 



write explicitly the six equations in (|3U|) as 

r-3-o- = XA1A2 (34) 

r^i-a = XA2A3 (35) 

r~3_^ ^ A^s^i (36) 

r^i-o- = XA4A1 (37) 

r-3-fT = XA4A2 (38) 

r^i-a = XAiAs, (39) 

and for convenience we consider separately the concave and convex 
cases. The concave case occurs when one particle is in the convex hull 
of the other three. In the convex case no particle is in the convex 
hull of the other three. Instead of formal definitions we appeal to two 
general examples of concave and convex conhgurations represented by 
our drawings in Fig. 1 and Fig. 2. 

A) First the concave case where a particle, labeled 4, is in the con- 
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vex hull of the other three, labeled 1, 2, 3. Assume with no loss of 
generality that 

Ai > ^2 > ^3 > > ^4 . (40) 

Then the first three equations (|34ll36p have the sign of A and the 
last three equations (j37ti39p have the opposite sign. Assume rij is 
the largest of (ri2, ^23, rsi), then the equilateral triangle based on r^j 
overlaps triangle 123. Therefore, since mass 4 is in the convex hull of 
the others, distances < r^j and hence r^f — a > r^^ — a. Since one 
of these is positive and the other negative, we can only have 

r^/ - a > , r^j' -a<^, (41) 

from which follows that 

A < . (42) 

As a Corollary of the previous proof, distances ryi , r23 , r^x are larger 
than o~'^l^ ^ which is larger than r4i,r42,r43. 

B) Second, the convex case where no particle is in the convex hull of 
the other three. We label the ends of the two diagonals with numbers 
1, 3 and 2, 4, respectively. The four labels are in the cyclic order 1234. 
Assume with no loss of generality 

^1 > A3 > > A2 > ^4 • (43) 

Thus equations (f36l) and ([38|) have the sign of A and the other four of 
the set (j34ti39p have the opposite sign. Assume Vij is the longest diago- 
nal (ri3 or r24) and rjfc is the shortest side of the set (ri2, r23, r43, r4i.) 
Then r^j > rj^ since assuming the contrary the four sides would be 
longer than r^j, and then the two equilateral triangles with com- 
mon side Tjj could be covered by the quadrangle 1234 and then r^j 
would not be the longest diagonal, contrary to the hypothesis. Then 
r^^ — a < r^^ — a , but since they must be of opposite signs we must 
have 

r-^-a<Q, rr,3-a>0, (44) 

from which 

A < . (45) 

As a corollary of this proof, the diagonals are larger than a~^f^ which 
is larger than the longest of the four sides. 

A different proof that A is negative is obtained in [3] as corollary 
of lemma 1. The sign of A is important because numerical solutions 
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to equation (j29p of the form (j30p may be found with A positive, that 
obey restriction (jlip but are not planar central configurations since 
they do not satisfy nor p^ . 

Additional relations between distances are obtained from the quo- 
tient, member by member, of equations (|34ll39p . One has 



(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
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— a 
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— a 
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— a 


'42 
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'23 


— a 


^^31 


— a 
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^23 


^31 
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^31 


— a 
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— a 


».-3 
^^31 


_^-3 
Ml 


A2 


'23 


— a 


^-3 
'42 


— a 


^-3 
'23 


'42 


A4 


Ml 


— a 


^-3 
M3 


— a 


'41 


'43 


A2 


M2 


— a 


'23 


— a 


^^12 


'23 




'42 


— a 


M3 


— a 


'42 


'43 


Ai 


M2 


— a 


'^31 


— a 


^-3 
^12 


^^31 



which are independent of the parameter A. The last expression on 
the right hand side of these equations results from the previous two, 
assuming that they are different. The middle terms imply restrictions 
that may be also obtained from the original equations, namely 

{rii - <T){rli -a) = {r^^ - a){r^i - a) = {r^i - a){rl^ - a) , (52) 

which all equal A^yliyl2^3^4- 

From these equations one also obtains the a parameter in terms of 
the distances 



r 3^ 3 
'12 '43 


'23 '41 


ly 3 3 3 3 
_ '31 '42 ' 12 '43 


^^12 + ^^43 


'23 '41 


3 1 3 3 3 
'31 '42 '12 '43 




T^~3 — 3 
'23 '41 


;v^~3 — 3 

'31 '42 




^2"3' + ^4"l' - 


'31 '42 



cr 



(53) 

Indeed this result is an identity after substituting the reciprocal cubed 
distances ([30|) . in the right hand side of this equation. 

All of equations (I46II53P appear in Schmidt [6]. We reproduce them 
here because they will be useful to prove further important results. 
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Coming back to the concave case we assume now the inequahties 
(cf. equation ([l0|) 

> ^2 > ^3 > > ^4 • (54) 
Then equations (j46ll51|) imply 

ri2 > rsi > r23 > a~^^^ > ns > > ni . (55) 
Returning to the convex case, now we assume 

^1 > yl3 > > As > ^4 ■ (56) 
Equations (j46H5ip then yield two inequalities for the four sides 

cr'^^^ > r23 > ri2 > Ui , a'^^^ > ras > > , (57) 

4 Symmetric configurations 

Important symmetric cases are also deduced from equations (|46ti5ip . 

The necessary and sufficient condition for a kite symmetry is the 
equality of two weighted areas A^s, for example iff 

Ai = A3 (58) 

then 

ri2 = r23 and r^i = r^s (59) 

and, furthermore, 

Si = S3 and mi = 7713 . (60) 

The kite symmetry imposes the above restrictions. The values of the 
other weighted directed areas A2 and A4 are arbitrary, but at least 
one of them must be of opposite sign to that of Ai. 

Imposing the previous symmetry twice leads us to two important 
symmetric cases. 

1) Equilateral triangle symmetry occurs iff 

Ai = A2 = A3, (61) 

which implies 

r23 = r3i = ri2 > = r42 = r43 . (62) 
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The three particles 123 are at the vertices of an equilateral triangle 
and particle 4 is at the center of this triangle, so that 

Si = S2 = S3 = -S4/3 , mi = 7712 = ma , — = ^/3■, (63) 

with A4 arbitrary, but satisfying the important equation 

Aa mi , ^ 

IT = -3— • 64 

2) Rhomboidal symmetry occurs iff 

^1 = A3 > > yl2 = ^4 , (65) 

implying that 

ri2 = r23 = r43 = r^i , and Si = S3 = -S2 = -S4 . (66) 

In order to yield the square symmetry, the diagonals must be equal, 
requiring the extra condition 

A1A3 = A2A4 . (67) 

3) Another kind of symmetry is obtained iff 

Ai = -A2 and A3 = -A4 (68) 
which imply the isosceles trapezium symmetry: 

r3i = r42 and r23 = (69) 

since 

A1A3 = A2A4 and A2A3 = A1A4 . (70) 
In such a case of trapezium isosceles symmetry we also find that 

Si = —S2 and S3 = —S4 , (71) 

which imply 

mi = m2 and ma = m4 . (72) 

Note that assuming the trapezium isosceles symmetry through (I69p . 
equation (j70p follows with the possible symmetries of the rhombus or 
of the isosceles trapezium. 

The isosceles trapezium symmetry gives the stronger square sym- 
metry if the four Aj^s have the same absolute value. 
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5 A new algorithm 



Equations ([TT]) . p3]) and ([3T|) have been often used by several authors 
to compute the six distances rji. from given values of the four masses. 
According to D. Saari [2] this perspective makes the problem difficult 
to manipulate. For example, with this perspective in mind, some ad- 
vances for the case of three equal masses with kite symmetry have 
been published by Bernat et al. [TO], with important results about bi- 
furcation values corresponding to different number of solutions. With 
our approach the problem with two equal masses is simpler than the 
Bernat et al. paper where the restriction of three equal masses leads 
to a more complex algorithm. 

Since the lengths and masses are defined up to arbitrary units and 
noting that the parameter a has units of the inverse of a volume, we 
assume with no loss of generality that all the distances are measured 
in the unit of distance a~^^^ which simplifies equation ([30]) into 

r7^3 = 1 + xAjAk . (73) 

In the following, we shall therefore use the simplification a = 1, and 
we also assume that the directed weighted areas Aj. are known as four 
given constants (the arbitrariness of the four values is justified by the 
arbitrariness of the four masses.) The previous equation gives all the 
distances as functions of the unknown parameter A and the four con- 
stants Ak- From them, using Heron's formula (llSp . the areas of the 
four triangles become functions of A. The sign of a directed area Sj 
is inherited from the same sign of the corresponding constant Aj. Re- 
striction ()14p then determines the value of A for a plane solution and 
from it, the value of the six distances and the four masses are found 
(in terms of the four constants Aj ) . This is an implicit way to deduce 
planar central configurations with four masses. In many non symmet- 
ric cases restriction (llip is sufficient to define a central configuration, 
but in a few cases, restriction (jl4p is not automatically satisfied and 
allows one to discriminate non-physical situations (such as negative 
distances or evidently non-planar solutions or geometrically impossi- 
ble figures.) We stress the fact that to determine A, we chose to find 
the root of (jl4p as a function of A. In any case considered numerically 
this solution is unique. 

This approach to the four body central configurations is similar 
(see [2]) to considering the Euler collinear central configurations of 
three particles, where instead of giving the ratio of the masses and 
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computing the corresponding ratio of two distances with a quintic 
equation, the ratio of two distances is first given and the ratio of the 
masses is computed from just a linear equation. 

The symmetric cases of kite symmetry, equilateral triangle sym- 
metry, rhomboidal symmetry, and trapezium isosceles symmetry were 
presented from this approach, based on the weighted directed areas 
A,. 

In general, given a solution, it is invariant to a change of sign of 
the four constants Aj. Therefore, with no loss of generality we may 
assume ^1,^43 positive and ^4 negative, henceforth the concave case 
is determined by A2 positive and the convex case by A2 negative. 

More generally, consider that a solution with a set of values of 
the four constants Aj is determined by the value of A. Multiplying 
the four constants by a real k: Aj — > kAj the solution requires a 
transformed A as: A — > X/k'^. 

This procedure produces in most cases a solution, but we find some 
cases where values of these constants yield none. Examples are given 
in the next section devoted to increase our study of the kite symmetry. 



6 Asymptotic cases of Kite Symmetry 

In order to deduce more of its properties, we come back to the kite 
symmetry recalling equations (j58ti60p that we rewrite here 

^1 = ^3 , 

then 

ri2 = r23 and = r^s 

and further 

= 53 and mi = . 

We shall consider four asymptotic cases of kite symmetry where 
one of the four masses approaches either zero or infinity. That mass 
lies on the symmetry axis between masses mi and m3. For clarity we 
subdivide these four cases in different subsections. 

We discuss first the concave solution with Ai = A3, A2 positive 
and A4 negative. m4 is in the convex hull of the other three particles. 
Particles with masses m2 and m4 are placed on the symmetry axis of 
particles with equal masses mi and m3. 
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Instead of the planar condition ()14p we use Pythagoras' Theorem 
which yields the equivalent equation 

-42-\/r?2-^ + \/r|i-^=0. (74) 

This constraint on the distances between particles is simpler than the 
other condition, resulting in a smaller numerical error, although it is 
only valid for the concave kite symmetry. 



6.1 Euler convex 

Numerical evidence was obtained considering fixed values of j4i and 
A4, and A2 variable. Starting with A2 = Ai = A3, the equilateral 
triangle solution, we decrease only the A2 value. We observed that as 
A2 is decreased, m2 approaches zero until an asymptotic value of A2 
is reached where there is no solution. Since A2 is a positive number, 
and the mass m2 approaches zero, this implies that the directed area 
52 also tends to zero, but this is possible only when r^i approaches 
the value r3i/2, which was verified numerically. Particles 1, 4, 3 be- 
come asymptotically collinear, with particle 4 at the midpoint of ri^. 
Imposing the equality rn = r3i/2 the allowed limiting value of A is 
found to be 

a function of Ai and A4 only. Substitution in the constraint (j74p 
as function of the constants Aj and A (that becomes the Pythagoras 
Theorem of the right triangle formed by 1, 2, 4,) one deduces the 
limiting minimum value of ^2/^1 as the root x of the equation 
1 11 

(76) 



2 



(8^1 - ^4 - 7^ix)3 (8A1 - ^4 - 7^42;)3 (8^1-8^4! 

This root is a function of the ratio —A4/A1 only. As this ratio ap- 
proaches zero, the limit of A2/A1 tends to 

^ (8 - \/8) = 0.738796125... . (77) 

On the other hand, as the ratio —A4/A1 approaches 00, the limit of 
A2/A1 tends to 

- ( -4= -1) = 0.077085817... . (78) 
7 VV27 J 
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In this limit, as the mass m2 tends to zero, the other three particles 
tend to lie on a straight line, a central configuration first studied by 
Euler. 



6.2 Lagrange concave 

In the second asymptotic case, we again start from the equilateral 
solution with Ai = = A2, all positive, and a fixed A4 < 0, but 
now we increase the value of A2- As A2 grows to infinity the three 
distances r4i,r43,r3i, approach 1, m2 tends to zero, ri2 = tend to 
a limit larger than 1, r42 has a limit smaller than 1, and A approaches 
zero from below. This limit for A implies that r4i,r43,r3i tend to 1. 
The peculiar behavior of ri2 = ^23 and r42 is explained by the property 
that 

lim A^2 = finite number < . (79) 

A— >-0,A2^oo 

To compute the unknown distances r42 and ri2 we may apply the 
first right hand side of equation (jSip with a = 1: 



A. 



4 



A, - 1 

In this limiting condition, the angle 143 is vr/S and angle 241 is Svr/G. 
The trigonometric cosine-sides relation for triangle 241 gives in the 
limit case the additional relation 

l + r|2 + r42\/3 = r?2- (80) 

Eliminating ri2 from these two equations we obtain an equation for 
r42 that is a function of the ratio A4/A1 



A4 _ rj - 1 



Ai (l + r22 + r42\/3) -3/2-1 



(81) 



In this C&SG clS A2 — > 00, the mass 1112 tends to zero. The other 
three particles are in an equilateral triangle central configuration first 
studied by Lagrange. Prom this value of r42, the limiting value Ayl2 
may be computed from equation (|38p as 
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6.3 Lagrange convex 

We now turn to the convex solution with Ai = A3 positive and A2, 
Ai negative, no particle is in the convex hull of the other three. Just 
as in the concave case, particles with masses m2 and are on the 
symmetry axis of the particles with equal masses mi and ms. 

The Theorem of Pythagoras in this case gives us the constraint 




r^^-fl2-'f-fli-'-f = ^. (82) 

valid only for the convex configuration with kite symmetry. 

We have studied numerically this convex case of kite symmetry 
keeping the two positive weighted areas Ai and A3 and the negative 
weighted area A2 constant and tabulating the six distances and four 
masses when different values of the weighted negative area are 
given. Any chosen negative value of this constant gives a solution. 

In the limit as tends to minus infinity, A goes to minus zero 
with the property 

lim AA4 = finite number > . (83) 

A^-0,A4^-oo 

In this limit, the mass 4 approaches zero and the other three particles 
form a Lagrange's equilateral solution with the three distances ri2, 
^"23, '"SI tending to the value 1. The unknown distances rn and r42 
may be calculated in similar form to the convex Lagrange case. The 
two distances are related by equation (f49|l 



M ^ r^t - 1 
A2 rli - 1 ■ 

Again a trigonometric relation between the two distances is obtained 
from the triangle 124. Angle 124 is 7r/6; the cosine-sides relation for 
this triangle produces 

l + r|2-r42\/3 = r|i. (84) 

Elimination of r4i between the last two equations gives side r42 as a 
function of the ratio A1/A2, as the root of the equation 

A, (l+r|,-r42^/3)-^/^-l 
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With this value of r42, the finite Umiting value of XA4 may be 
computed now from equation (|38p 



AA4 



r 



-3 

42 



- 1 



A2 



6.4 Coorbital 



Turning presently to the examination of the convex solution with Ai 
= A3 positive and A2, A4 negative in the limit where the negative 
constant A4 tends to zero, we find that the three distances r4i, r42, r43 
approach the same limit 1. In this limit, the three particles 1, 2, 3 
lie on a circle with center at the position of particle 4, with a mass 
tending to infinity. This corresponds physically to the coorbital case of 
particles of negligible mass, around a central large mass, whose study 
was pioneered by Maxwell in 1856. See for example reference [H] for 
a panoramic view of this problem. 

The kite symmetry implies 9 = angle 142 and that 29 =angle 143. 
Therefore the two distances ri2 and ri3 are given by 



They are related by the trigonometric identity 



Distances ri2 and r^i satisfy also the equation ([T7|) with a = 1 



Eliminating r^i between the two last equations allows us to deter- 
mine ri2 from the equation 




(87) 



A2 _ r^i - 1 
As - 1 



A2 

A3 




(88) 



The limiting value of A, as IA4I — > 0, is deduced from 



A = 
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In comparison, the approach used in [14j is to assume the restric- 
tion that the three satehite masses are equal. Then two non trivial 
solutions are given by those authors in terms of angle 9, one corre- 
sponding to the convex case, and the other to the concave case. Ac- 
tually, we find that these angles should both obey the trigonometric 
restriction 



1 



1 



|sin'^(0/2)| 



+ 2 cos (6*) 



1 



1 



8|sin^(0)| 



0. 



with the two non trivial solutions 



e = 0.826602936080376.. , 6 = 2.4219145305912.. . 

To relate with our approach we stress that we start with the values 
for the weighted directed areas 



= ^3 = 1 , ^4=0 , A2 



12 



1 



(2sin(6'/2) 



-3 



1 



'31 



(2 sin ( 



1 



Then, computing the masses we found that masses mi=m3=m2, as 
might have been expected. Thus the relation of our method with [14] 
is established for the (1 + 3) central configuration. 



7 Four different masses 

In the previous section, several limiting cases were presented where 
some symmetric limit central configurations were computed from sim- 
ple expressions and compared with numerical experiments. 

In the general case of four different masses we have not much the- 
ory to add to the previous existing Dziobek's [4] findings. We have 
nevertheless to offer the practical algorithm presented in this paper 
where the four constants Aj are assumed known and from them the 
distances, and hence the areas, are known as functions of A that are 
determined from a plane constraint. At the next stage, the four masses 
are determined instead of giving them from the beginning. 

The case with four different masses and no evident symmetry has 
also asymptotic limiting configurations similar to those cases described 
in the previous section supported by symmetry and is dealt with using 
the same new algorithm. 

First the case where one of the weighted directed areas, say 
tends to infinity, and simultaneously the A parameter goes to zero. 
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The three distances that are expressed in terms of three finite, non- 
zero, values of the directed weighted areas approach unity, because of 
the A behavior. They may be assumed 1 from the beginning. Intro- 
ducing a new finite parameter equal to the formal product XA^, this 
is now determined from the plane constraint. This limit corresponds 
to the physical situation in which one mass disappears and the other 
three become a Lagrange's equilateral three body central configura- 
tion, which is completely determined by the equilateral condition. 

The other central configuration of three bodies, the Euler collinear 
configuration, occurs again in the limit where the fourth mass tends to 
zero. In this case both the mass and the area associated to the vanish- 
ing particle approach zero; the quotient, the corresponding weighted 
directed area takes a finite non-zero value. To compute the limiting 
case we assume that only three weighted directed areas of the three 
collinear particles are known, then introduce the collinear condition 
on the three distances, function of A, between the collinear particles 
in order to obtain the A value. That this is the usual Euler condition 
comes from Albouy's proof [3] of Dziobek's equation ([30]) that is also 
valid for the collinear three body case with the directed weighted areas 
replaced by directed weighted distances. Because the three triangles 
have the same height, directed areas and directed distances are pro- 
portional. The unknown value of the fourth weighted directed area 
is determined from the planar constraint ()14p because it is the only 
remaining unknown in the expressions for the six distances. This lim- 
iting central configuration can be obtained from a concave solution 
when one of the weighted directed areas changes its value. The de- 
duced limiting value of ^4 bounds this weighted directed area from 
further change. 

The coorbital case when one of the masses tends to infinity is 
obtained in the zero limit of one particular weighted directed area, 
assuming the other three weighted directed distances are known, finite 
and non zero. Then the three distances connecting the massive particle 
become equal to one. We may replace this value from the beginning. 
The other three weighted directed areas are known in terms of A, which 
is determined as usual. These three distances determine the geometric 
angles separating the three satellites (with three different masses.) 

In the generic case of different masses, the four assumed known 
weighted directed areas have been related to one of the coordinates 
of the four particles determined by the internal rotation around the 
center of mass of a rigid tetrahedron, which is a function of the four 
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masses only, having equal principal moments of inertia. This is the 
kernel of the coordinates introduced by Piha to describe the general 
dynamics of the Four-Body Problem [15] that will be presented in a 
separate paper. As this kernel rotates (by three Euler angles) it then 
collapses into a plane, loosing the four coordinates in the direction or- 
thogonal to the plane. These four coordinates, function of the masses 
and two Euler angles, form four quantities that are proportional to 
the four weighted directed areas. We refer to that paper for details 
about these coordinates. 
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